EPJ manuscript No. 

(will be inserted by the editor) 



OO 
O 

o 

(N 

> 

o 
O 



> 
OS 

m 

(N 

d 

00 

o 

> 

>< 



Power laws in elementary and heavy-ion collisions 
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Abstract. We review from the point of view of nonextensive statistics the ubiquitous presence in elementary 
and heavy-ion collisions of power-law distributions. Special emphasis is placed on the conjecture that 
this is just a reflection of some intrinsic fluctuations existing in the hadronic systems considered. These 
systems summarily described by a single parameter q playing the role of a nonextensivity measure in the 
nonextensive statistical models based on Tsallis entropy. 



PACS. 89.75.-k complex systems - 24.60.-k statistical theory and fluctuations 
duction (relativistic collisions) - 25.75.-q relativistic heavy ions collision 



25.75.Dw particle pro- 



1 Introduction 

In many domains of physics, especially in elementary and 
heavy-ion collisions, for decades the prevailing understand- 
ing was that exponential shape of most of the observed 
spectra of produced secondaries suggests their statistical 
(or even thermodynamical) origin. Therefore, when look- 
ing for spectra of transverse momenta it is commonly 
assumed that the observed inverse mean transverse mo- 
menta (characterizing the widths of such exponential spec- 
tra) play the role of a temperature of the hadronizing sys- 
tem [I]. This assumption allows us to use, in what fol- 
lows, the whole machinery of statistical models. This is 
specially important in the case of high energy collisions of 
heavy ions because it allows us to use the tools of statisti- 
cal physics to investigate the possibility of the formation 
in such collisions of a new state of matter, the so called 
Quark Gluon Plasma (QGP) [2] (in this case one is in- 
terested in details of hadron -<=>• QGP phase transition, 
which can only be investigated in this way). 

However, one should always keep in mind another pos- 
sibility, namely that such behavior can just be due to the 
fact that, when presenting our data, out of many particles 
produced (tens or thousands at present, say N to be spe- 
cific) we select only one to make the corresponding plots. 
Single particle distributions averaged over many events 
are what is usually published. But this means that the re- 
maining (N — 1) particles will act as a kind of heath bath. 
Assuming that this heath bath is homogenous and large it 
is natural to expect that its action can be described by a 
single parameter, which we then call temperature, T, and 



identify with the temperature encountered in statistical 
models [3]. What one apparently observes then is just the 
usual Boltzmann-Gibbs (BG) statistics at work. 

The above reasoning assumes that any dynamics of the 
set of remaining (N — 1) particles is mostly averaged and 
what results looks very much like a state of hadronic mat- 
ter remaining in thermal equilibrium characterized by a 
temperature T. Some effects, however, can survive this 
equilibration process and can show up as apparent de- 
partures from the assumed thermal equilibrium. This is 
usually regarded as a departure from BG statistics and, 
finally, considered as a kind of failure of the simple statisti- 
cal approach. Such observations are therefore a subject of 
separate investigations in which many different (dynam- 
ical) ideas compete (like, for example, the production of 
resonances and the flow of matter, to name only two) [4]. 
However, because there is only a limited amount of avail- 
able data, usually one cannot decide which of the proposed 
dynamical remedies is right or, if we agree that they should 
all be present, in what proportion they show upp 

On the other hand, one can argue that, perhaps, it is 
the form of statistical model used which should be modi- 
fied in such a way as to account (at least to some extent) 
for detected irregularities (i.e., for departures from the 
BG approach). Therefore, instead of inventing and inves- 
tigating different dynamical assumptions, one can instead 
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1 This is connected with the important problem of how much 
information given measurements are providing us with. This 
is, so far, only sporadically discussed in the domain of high 
energy multiparticle production processes using information 
theory approach in which the entropies mentioned here are re- 
garded as measures of information [5][6]. This subject is, how- 
ever, beyond the scope of our review. 
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investigate the possibility of replacing the usual statistical 
model based on BG entropy, Sbg, by its modified version 
based on some other form of entropy discussed in the lit- 
erature [7] . Such models are widely known nowadays from 
other branches of physics and are used whenever a phys- 
ical system under investigation shows memory effects of 
any kind, experiences long range correlations (i.e., is in a 
sense " small" because its size is comparable with the range 
of forces acting in it), experiences some intrinsic fluctua- 
tions, or the phase space in which it operates is limited or 
has fractal structure. In this review we shall discuss the 
application of such a model taken in the form proposed 
by Tsallis, i.e., in the form based on the so called Tsallis 
entropy, St = S q [8]: 



(called also Tsallis distributions 



S a 



9-1 



Sbg 



Pi In. pi 



= i. The S q is nonextensive because 



Notice that Sbg = S q 

for any two independent systems A and B (in the usual 
sense, i.e., for which Pij(A + B) = pi(A)pj(B)), one ob- 
serves that 

S q (A + B)= S q (A) + S q (B) + (1 - q)S q (A)S q (B), (2) 

In this sense the entropic index q is a measure of the 
nonextensivity in the system without, however, directly 
showing up its cause. This must be provided from else- 
where. 

To shed more light on the physics involved here, let 
us come back to the previous reasoning with some effec- 
tive thermal bath being formed by the (N — 1) particles 
remaining after selection of the one used for making final 
histograms. Notice that, in high energy multiparticle pro- 
duction reactions observed in elementary and heavy-ion 
collisions we are interested in here, such thermal baths 
usually (i.e., after more detailed scrutiny) do not satisfy 
conditions allowing us to introduce the notion of ther- 
mal equilibrium in the BG sense: they are always finite 
and can be hardly considered as being homogenous (in 
fact, in many cases they occupy only a fraction of the al- 
lowed phase space [5] or even have a fractal-like structure 
|10j and it is known that usually the hadronizing system 
under consideration experiences long range correlations). 
This means therefore that such a heath bath cannot be de- 
scribed by a single parameter T. The simplest thing to do 
seems to be to allow for some fluctuations of the parameter 
T and to replace it by its mean value T — > To = (T) and 
by one more parameter describing its fluctuations, using, 
for example, the (normalized) variance [TT1I12) : 



<(*: 



'J_\2 



(3) 



In this approach 



q = 1 + lu. (4) 

It can be next shown that such a heath bath leads in 
a natural way to the following ^-exponential distribution 



exp 



X 

T 



exp 



X 

T 



i-(i-«) 



(5) 



This is the power law we were searching for in differ- 
ent reactions [H V T5 U T5 V TH V T7 y TS V Tg y ^n | I^T V ^S] . the results 
of which will be reviewed herqj 

The physical picture presented above can be made 
more formal by saying that one replaces here the notion 
of strict local thermal equilibrium, customarily assumed 
in all applications of statistical models, by the notion of 
some kind of stationary state, which is being formed in 
the collision and which already includes some interactions. 
This concept can be introduced in different ways. For ex- 
ample, in |29j it was a random distortion of energy and 
momentum conservation caused by the surrounding sys- 
tem which resulted in the emergence of some nonextensive 
equilibrium. In [30U31] the two-body energy composition 
in transport theory formulation of the collision process is 
replaced by a generalized energy sum, h(E\,E2), which 
is assumed to be associative but which is not necessarily 
simple addition and contains contributions stemming from 
pair interaction (in the simplest case). It turns out that 
under quite general assumptions about the function h, a 
division of the total energy among free particles is possi- 
ble. Different forms of the function h then lead to different 
forms of entropy formula, among which one encounters the 
known Tsallis form. The origin of this kind of thinking can 
be traced back to the analysis of the g-Hagcdorn model 
proposed some time ago in (321. 

We close this section with a historical note. The recog- 
nition that some, apparently unexpected power law dis- 
tributions can be due to fluctuations came to us from the 
observation in cosmic ray physics [33j that there exists a 
long flying component (LFC) phenomenon in the propaga- 
tion of the initial flux of incoming nucleons. For example, 
instead of the normally expected exponential fall off of the 
depth distribution of the starting points of cascades, z, 



dN(z) 
dz 



const • exp (— z/X) . 



(6) 



one rather observes an exp g (— z/X) distribution, i.e., Tsal- 
lis like power law behavior given by Eq. ([5]) with z re- 
placing X (here A ~ 1/a is the mean free path describ- 
ing propagation of the incoming flux in the atmosphere 
with a being the relevant cross section). On the other 
hand, some time ago we have shown in [33] that this effect 



2 Except for [T2] discussion concerning the meaning of the 
parameter q was limited to the case of q > 1 only. As already 
mentioned in [12] , the case of q < 1 seems not so much con- 
nected with any genuine fluctuations but rather, in some way, 
with limitations of the allowed phase space [13] . It is worth 
mention that the idea of possible fluctuations of otherwise in- 
tensive quantities has already been formalized by introducing 
a new concept of so called superstatistics [14] . 

3 It should be noticed that there are also some other investi- 
gations on this subject |~~ 



] , with rather similar 
conclusions but which we will not discuss here. 
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can be explained by assuming that hadronic cross sections 
should be regarded as fluctuating quantities with widths 
(defined as normalized dispersion), w a = (a 2 ) /(a) 2 — 1 
(and growing logaritmically with energy) . We were at that 
time prompted by the fact that such an idea was widely 
investigated in the usual hadronic collisions [35J3- It was 
then quite natural to connect the parameter q with fluc- 
tuations. This done in [111112] , as mentioned above. 

We shall review our results in this field in the next 
Section. Section [3] contains some new recent developments 
in this field. The final Section contains our conclusions and 
a summary. 



2 Review of fluctuations in multiparticle 
production processes 

High energy collisions result in a multitude of particles of 
different kinds being produced. Most are just mesons of all 
kinds (overwhelmingly pions). For those who are looking 
for some new and/or rare phenomena they form unwanted 
background which must somehow be substracted, for oth- 
ers they are a subject of thorough investigations allowing 
us to look inside the very early stages of the collision pro- 
cess as well as at the hadronization stage of the matter 
produced (proceeding probably via the formation of the 
QGP, for example). In both cases a simple and trustwor- 
thy representation of data is very important, this justifies 
our investigations in this field to be reported here. 



2.1 Generalized heat bath 
temperature 



fluctuations of 



We first recall the physical picture behind the generalized 
heat bath introduced in Section [T] which we have proposed 
in [TT|ll2j. Our reasoning was as follows. Suppose we have 
a thermodynamic system, in a small (mentally separated) 
parts of which the temperature can take different values, 
i.e., in the whole system it fluctuates with AT ~ T. Let 
£ (t) describes stochastic changes of temperature in time. If 
the mean temperature of the system temperature is (T) = 
Tq then, as a result of fluctuations in some small selected 
region, the actual temperature T" equals 



T' = T - b€(t)T, 



(7) 



where the constant b is defined by the actual definition 
of the stochastic process under consideration, i.e., by £(£), 
which is assumed to satisfy the condition that 



(£(*)) = o 



(8) 



and which correlator, (£(t) £(t + At)), for sufficiently fast 
changes is equal to 



(£(*)£(* + At)) = 2D5(At). 



(9) 



4 It is interesting to notice that this idea has been revitalized 
very recently in [36] and connected with the fluctuations in the 
gluonic content of hadrons. 



The inevitable exchange of heat between any selected re- 
gion of the system and the rest leads to equilibration of 
the temperature in the whole system. The corresponding 
process of heat conductance is described by the following 
equation [37], 



dT 

c pP — - a(T' - T) = 0, 



(10) 



where c p , p and a are, respectively, the specific heat under 
constant pressure, density and the coefficient of external 
conductance. Using T' as defined in ([JJ we finally get the 
linear differential equation for the temperature T with r 
b- 



dT 
~dt 



•£(*) 



T 



1. 



-T . 



(11) 



It can be now shown that this equation leads to the Lange- 
vin equation with multiplicative noise term resulting in 
fluctuations of the temperature T given in the form of a 
gamma function [11] 



f(T) 



and characterized by the parameters p and a, 



r(a) 



M = 



D 



and 



1 

7d' 



(12) 



(13) 



This is to be compared with Eq. (J3J in which now lo = tD. 
Function f(T) as given by Eq. (fT2"]) is the distribution that 
should be used to smear the parameter T in the usual 
exponential distribution of the BG statistical model and 
which results in the Tsallis distribution, Eq. (0. To sum- 
marize: a small addition of the multiplicative noise de- 
scribed by a damping constant in the Langevin equation 
results in a stationary distribution of particle momenta, 
which develops a power-law tail at high valued- 



2.2 Transverse and longitudinal dynamics 

To begin our presentation we first set the stage. The char- 
acteristic pattern of the multiparticle production processes 
is that most of the secondaries are produced with small 
transverse momenta pr (mostly below 1 GeV) and are 
therefore concentrated in the longitudinal phase space given 
by the longitudinal momenta ph (which is described in 

+ Pl 



terms of the rapidity y = ^ In gjr^ , where E = 



with m,T = y m 2 + p\, being the so called transverse mass 
(m is mass of the particle); in other notation E — mx cosh y 
and pl — rriT smhy). The terms transverse and longitu- 
dinal are defined with respect to the direction of the col- 
liding particles. Data are presented as distributions either 



Actually, as shown in |14] when discussing super statistics, 
there is a whole class of functions leading from exp(X) to 
exp q (X). But only this has a simple physical interpretation 
as presented here. More general version of Langevin equation 
containing also additive noise have been considered in [38] . 
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in px or in y. In both cases they show exponential behav- 
ior either in pt or in the energy E — (%) coshy (with 
yjm? + (p T ) 2 ): 



{m T ) 

dN 
dp T 



= C Pt pt exp ( - 



Pt 
T 



dN 
dy 



= C y exp - 



(14) 



One observes dramatic differences in both distributions 
reflected by the differences in the values of the parameter 
T, which is of the order of one hundred MeV in px space 
(where T = T PT and is universal, i.e., essentially energy 
independent) and tens of GeV (depending on the energy of 
collision) in p^ (or y) space (where T = T PL and depends 
on energy). This means that the two distributions reflect 
different physics: those in pt space are believed to be es- 
sentially " thermal- like" and subject to a thermodynamic 
interpretation whereas, those in p^ space are sensitive to 
the available energy and to the multiplicity of produced 
secondaries. Because of this their fluctuation patterns will 
be different, i.e., when described by Tsallis power-like form 
Eq. §5§ the corresponding parameters (qr — 1) and (qz — 1) 
will differ dramatically. Also the physical meaning of these 
parameters will be different reflecting different sources of 
fluctuations. 



2.3 Longitudinal phase space 

We start with the longitudinal distribution in rapidity (av- 
eraged over px)- In Fig. [I] one observes that q < 1 is ef- 
fectively cutting off the allowed longitudinal phase space 
(here defined by the initial available energy M = 100 GeV 
and assumed constant transverse mass wit = 0.44 GeV 
and weakly depending on the assumed multiplicity of the 
produced particles N) . Actually, from Eq. ([5]) it is obvious 
that only such combinations of q and X and A are allowed 
for which [1— (1 — q)X/X] > 0. In [13], when fitting longitu- 
dinal distributions without restricting the available energy 
by introducing the so called inelasticity coefficient K < 1, 
the only role of q, which was found to be q < 1 there, was 
to limit the amount of energy used (showing the necessity 
of introducing inelasticity when considering multiparticle 
production processes, cf., [33] for review on this subject). 
For q > 1 one observes a visible enhancement of distribu- 
tion tails. 

Physica A 344, 568 (2004) 

The examples of fits to the actually observed single 
particle distributions in rapidity are shown in Figs. [2] and 
0] In the left panel of Fig. [2] (see [TH] for details) results for 
pp and pp collisions at energies varying between y/s = 20 
GeV to 1800 GeV are displayed. From each listed energy 
of collision, E cm = y/s, only a fraction K q has been used 
for the production of secondaries (according to [33] )■ The 
other input was the mean multiplicity of charged secon- 
daries produced in nonsingle diffractive reactions at given 
energy: fi c h = — 7.0 + 7.2s 0127 fTJ (corresponding to the 
total number of produced particles, N — %n c h)- The al- 
lowed phase space is one dimensional with only a small 
energy dependence of the mean transverse momentum al- 
lowed, (p T ) = 0.3 + 0.044 In (Vs/20) PQ (all secondaries 
will be assumed to be pions of mass p = 0.14 GeV). 




Fig. 1. Examples of the most probable rapidity distributions 
as given by Eq. (|14|l for hadronizing mass M — 100 GeV de- 
caying into N secondaries of (transverse) mass mr = 0.4 GeV 
each for different values of parameter q = Ql (reproduced by 
permission of Springer- Verlag from \18\). 

As discussed in detail in [19] , one gets in this case not 
only the parameter q but also the true inelasticity, K, of 
the reaction (in fact, even, for the first time, its distribu- 
tion, x(K)). Let us, however, concentrate on the parame- 
ter q, which bears information on fluctuations. It turns out 
to be energy dependent as presented in Fig. [3] Surprisingly 
enough it turned out that the same behavior is observed 
for the inverse of k characterizing the so called Negative 
Binomial distribution (NB) 1 of the multiplicity of ob- 
served secondaries, which depends on two parameters: the 
mean multiplicity (n c h) an d the parameter k (k > 1) af- 
fecting its width (a)n c h is dispersion), 



I _ a 2 (n ch ) 
k (n ch ) 2 



1 



{rich) 



(15) 



For k — > 1 NB approaches a geometrical distribution where- 
as for A; - 1 — > it approaches a Poissonian distribution. In 
general it is found PQ that ± = -0.104 + 0.058 • ln^/s, 
which fits the obtained values very nicely, cf. Fig. [31 

To fully understand the possible physical meaning of 
the parameter q(— q^) in this case let us remind our- 
selves that, in general, the nonextensivity parameter q 
summarizes the action of several factors, each of which 
leads to a deviation from the simple form of the exten- 
sive BG statistics, as was mentioned before, out of which 
we are interested most in the possible intrinsic fluctua- 
tions existing in the hadronizing system [11] , Notice that 
in our fits we have not explicitly accounted for the fact 
that each event has its own multiplicity, N, but we have 
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a E m = 53 GeV 

* E™ = 200 GeV 

* E°[ = 540 GeV 
o E m - 900 GeV 
V E™ = 1800 GeV 

* E_ = 20 GeV 





Fig. 2. Examples of applying the nonextensive approach to 
longitudinal distributions. Left panel: fit to rapidity spectra 
for charged pions produced in pp and pp collisions at differ- 
ent energies [40]. Right panel: rapidity spectra measured in 
e + e~ annihilations at 91.2 GeV [41] (dotted line is for K q — 1 
and q — 1 whereas the full line is our fit with K q — 1 
and q — 0.6). (Reprinted from Physica A344, F.S. Navarra, 
O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, "Information the- 
ory in high-energy physics (extensive and nonextensive ap- 
proach)", 568, Copyright (2004), with permission from Else- 
vier; http://www.elsevier.com ). 



0.35 




» q-1 
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10 10' 
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Fig. 3. The values of the nonextensivity parameter q obtained 
in fits shown in the left panel of Fig. [2] compared with the 
values of the parameter k of a Negative Binomial distribution 
fit to the corresponding multiplicity distributions as given in 
[T|. ( Reprinted Fig. 6 with permission from F.S. Navarra, O.V. 
Utyuzh, G. Wilk and Z. Wlodarczyk, Phys. Rev. D 67, 114002 
(200 3). Copyright (2003) by the American Physical Society; 
URL |http://link.aps.org/abstract/PRD/v67/ell4002[ DOI: 
10.1103/ PhysRevD. 67.114002.). 

used only its mean value, (iV), as given by experiment 
where (N) = Y,N P(N) with P(N) being the multiplic- 
ity distribution Q On the other hand, the parameter T in 
this case is not so much a temperature, but only a kind of 
"partition temperature", understood as mean energy per 
produced particle, i.e., T ~ W/(N) (where W = Ky/s, 
where -^s is the total energy of collision) [33] • Therefore 



6 Actually, we have used only its charged part, (n c h)i assum- 
ing that N = |(n c h), i.e., neglecting in addition also possible 
fluctuations between the number of charged and neutral sec- 
ondaries. 



in this case one can just as well speak about the fluc- 
tuations of (N). Following therefore the ideas of [TT] we 
would like to draw attention to the fact that the value of 
may also be understood as the measure of fluctua- 
tions of the mean multiplicity (for example, in the usual 
Poissonian multiplicity distribution characterized just by 
a single parameter, the constant mean multiplicity n), and 
in the case when such fluctuations are given by a gamma 
distribution with normalized variance D(h), one obtains 
the Negative Binomial multiplicity distribution with 



I = D W = -four- 



This is because in this case one has 



(16) 



P(n) 



dn 



r{k- 



e n n n r y k n k 1 e 



n! 

n) 



r(k) 



r(i + n )r(k) ( 7 + i) fe +»' 



(17) 



where 7 = ^|y. 

Therefore the situation in longitudinal phase space is 
following: When there are only statistical fluctuations in 
the hadronizing system one expects a Poissonian form 
of the corresponding multiplicity distributions. The exis- 
tence of intrinsic (dynamical) fluctuations means that one 
allows the mean multiplicity fi to fluctuate. It is natural to 
assume that these fluctuations contribute predominantly 
to the longitudinal phase space, i.e., that D(n) = q — 1 
and that 

(18) 

This is observed in the data. 

The right hand panel of Fig.[2]displays results for e + e~ 
annihilations for which, by definition, K q — 1 (because al- 
ways all the energy of initial leptons is available for the 
production of secondaries) and which can be fitted only 
with q < 1 (in our case q = 0.6). This should be contrasted 
with results obtained describing the pt distributions in- 
stead where one finds q > 1 [23] . This point deserves closer 
scrutiny. The result for q = 1 clearly shows that observed 
discrepancies are not connected with the particular value 
of q, but rather with some additional mechanisms oper- 
ating here, the action of which would, however, change 
our results only slightly (for example, a possibility of two 
rather than one source or y-dependent (pr), as mentioned 
already in [6]). With the above reservations, let us then 
take a closer look at the possible origin of q < 1. We 
have already encountered a similar situation when in [13] 
q < 1 was simply closing the allowed a priori phase space, 
acting therefore as inelasticity parameter K. When con- 
sidered as a signal of fluctuations (similar to the q > 1 
case) [12] it causes trouble because in this case the tem- 
perature T does not reach an equilibrium state, in fact 
one now has that the source term (right hand side of Eq. 
([TT]) ) is T /r - (q- 1)E/t rather than T /r used for for 
q > 1 case (cf., [12] )• This means than that in this case 
we have a kind of dissipative transfer of energy from the 
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region where (due to fluctuations) the temperature T is 
higher. It could be any kind of convection- type flow of 
energy; for example, it could be connected with the emis- 
sion of particles from this region (for example, in our case 
from a quark (q) and antiquark (q) jets formed in the first 
e + e~ — * q + q to gluons and qq pairs and later on to fi- 
nally observed hadrons). This means that q < 1 signals 
that in the reaction considered, where K q — 1 and where 
we have to account for the whole energy exactly, conserva- 
tion laws start to be important and there is no possibility 
for a stationary state with constant final temperature to 
develop. Instead, the temperature T depends on the en- 
ergjQ, and for large energies tends to zero (notice that in 
this case one has a limitation on the allowed energy of the 
produced secondaries: E < T /(l — q)). This is not the 
case for the pr distribution analysis [53] because most pr 
are small and are not influenced by the conservation laws 
but instead reflect a kind of stationary state with q > 1. 



PHOBOS 
0-6% 




Fig. 4. Examples of applying a nonextensive approach to 
longitudinal distributions. Left panel: fits to NA49 data for 
Pb+Pb collisions [42]. Right panel: fits to PHOBOS data for 
Au+Au collisions [43] . (Reprinted from Physica A340, F.S. 
Navarra, O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, "Informa- 
tion theory approach (extensive and nonextensive) to high- 
energy multiparticle production processes", 467, Copyright 
(2004), left panel, and from Physica A344, F.S. Navarra, O.V. 
Utyuzh, G. Wilk and Z. Wlodarczyk, "Information theory in 
high-energy physics (extensive and nonextensive approach)", 
568, Copyright (2004), right panel; both with permission from 
Elsevier; http://www.elsevier.com ). 

Now look at the left hand panel of Fig. 2J It shows 
fits to NA49 data [42] on n~ production in PbPb colli- 
sions at three different energies per nucleon. The obtained 
values of nonextensivities and the corresponding inelas- 
ticities, (q;K q = 3 • Kf) are: (1.2; 0.33) for 17.3 GeV, 
(1.164; 0.3) for 12.3 GeV and (1.04; 0.22) for 8.6 GeV. The 
origin of q > 1 in this case is not yet clear. The inelastic- 
ity seems to grow with energy. It is also obvious that, for 
higher energies, some new mechanism starts to operate 
because we cannot obtain agreement with data using only 



7 Actually, in the case considered in |12j fluctuations depend 
on energy in the same way leaving the relative variance uj con- 
stant and leading to q = 1 — to. 



energy conservation. The best fit for 17.3 GeV for NA49 
data actually for the case of q = 1 and two sources sepa- 
rated in rapidity by Ay = 0.83 (cf., [6] for other details). 

Finally, the right hand panel of Fig. [4] presents fits 
to pion production in Au+Au collisions 43J for the most 
central events (covering collisions proceeding with impact 
parameter range — 6%) 0. They can be fitted by choosing 
K q = 1 and then q = 1.29, 1.26 and 1.27 for energies 19.6, 
130 and 200 GeV, respectively (cf., [20] for other details). 
As before, the origin of q > 1 in this case is not yet clear. 

Although the situation in AA collisions is not yet clear, 
we are quite confident that interpretation of the q param- 
eter offered here remains valid. But, before settling this, 
one point has to be addressed. Namely, the above q were in 
fact qL responsible for the longitudinal dynamics only. On 
the other hand, multiplicity distributions are sensitive to 

P = vPl + Pt an d' as we nave seen here, both ph and pr 
show traces of fluctuations by leading to q > 1. However, 
as we shall see below, (q T — l) « (q^ — l) (what fits nicely 
the fact that pr space is very limited in comparison to pl 
one). Because there are no data measuring px distribu- 
tions at all values of rapidity y, i.e., providing correlations 
between parameters (T; q) = (Tl; q^) for longitudinal mo- 
menta (rapidity) distributions and (T;q) — (Tt;<?t) for 
transverse momenta distributions, we offer only the follow- 
ing approximate answer. Noticing that q— 1 = a 2 (T)/T 2 
(i.e., it is given by fluctuations of total temperature T) 
and assuming that u 2 {T) = ct 2 (Tl) + <t 2 (Tt), one can es- 
timate that the resulting values of q should not be too 
different from 



q = 



q L T 2 + q T Tr 

rp2 



n 



T 2 



rp2 



1, 



(19) 



which, for Tl ^S> Tt, as is in our case, leads to the result 
that q ~ <7l, i.e., it is given by the longitudinal (rapidity) 
distributions only. 



2.4 Transverse phase space 

As discussed before, transverse phase space seems to be 
mainly dominated by the thermodynamical-like effects gov- 
erned by the temperature T [T]. It is therefore the best 
place too look for any fluctuations of temperature, i.e., to 
look for any deviation of the the inverse slope of trans- 
verse momenta distributions, dN/dpx, from an exponen- 
tial shape. That such deviations are really observed is seen 
in Fig. On the left hand panel we can see fits to pr spec- 
tra measured by the UAl experiment |47j in pp at different 
energies using Tsallis distribution, Eq. ([5]), with X/X — » 
Pt/T and with the following values of (T = It [GeV], q = 



Actually these data are originally presented not for the 
rapidity y defined by the energy E and the the longitudinal 
momentum pl but for the so called pseudorapidity r\ defined 
by the total momentum p and the longitudinal momentum pl 
instead. There is therefore some ambiguity when transferring 
them from r] toy because of the pour knowledge of the rapidity 
dependence of the mean transverse momentum needed for such 
operation. 
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q T ): (0.134, 1.095), (0.135, 1.105) and (0.14, 1.11) for 
energies 200, 500 and 900 GeV, respectively (the values 
of the parameter q obtained in analysis of transverse mo- 
menta in elementary e + e~ reaction is similar [23j ) . These 
values should be compared with the corresponding values 
of (T = Tl; q — q£) previously observed for rapidity dis- 
tributions, which are equal to, respectively: (11.74 , 1 .2), 
(20.39, 1.26) and (30.79, 1.29) at comparable energies, cf. 

ma. 



the system [48] 



E dV/dp 1 


UA1 




900 GeV 










200 Gov" 






540 GeV 


) 1 2 3 


4 5 6 7 8 9 




S+S->h+X p lat , = 200A GeV 
3.0 <y <4JJ 
■ NA35 

1.0 T= 191 MeV 

— o— q= 1.015T = 1&6 MeV 



P T [GeV] 



o.e OS 



P T [GeV] 



Fig. 5. Examples of applying a nonextensive approach to 
transverse momenta distributions. Left panel: fits to pr spec- 
tra from pp UA1 experiment [46] for different energies (see 
text for details (reprinted from Physica A340, F.S. Navarra, 
O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, "Information the- 
ory approach (extensive and nonextensive) to high-energy mul- 
tiparticle production processes" , 467, Copyright (2004) , repro- 
duced by permission of Elsevier, http://www.elsevier.com ). 
Right panel: fits to 5 + S data from [47] (reproduced by per- 
mission of IOP Publishing Ltd from [15|h 

The right hand panel of Fig. [5] shows an example of 
similar behavior observed for nuclear collisions. Such col- 
lisions are of special interest as they are the only place 
where a new state of matter, the Quark Gluon Plasma, 
can be produced [2] and, because of this, they are in- 
tensively investigated using a nonextensive approach (see, 
for example, [i"2 ] l2i p 24 p 25 | l28]). As one can see, the best 
fit is obtained for q > 1, albeit in this case the value of 
(q — 1), which is the real measure of fluctuations, is no- 
ticeably smaller than in the case of elementary reactions 
mentioned above. On the other hand, although very small 
(\q — 1| ~ 0.015), this deviation leads to a quite substan- 
tial relative fluctuations of the temperature existing in the 
nuclear collisions, namely one gets that AT jT ~ 0.12. 

The question then arises: if this is treated seriously, 
what we are really measuring, what physical observable 
does it correspond to? It is important to stress that these 
are fluctuations existing in small parts of a hadronic sys- 
tem in respect to the whole system rather than of the 
cvcnt-by-event type for which, 

AT/T = 0.06/ViV -> 

for large N. The answer is that the measured fluctuations 
provide a direct measure of the total heat capacity C of 



1 

c 



(20) 



(/? = i) in terms of w = q — 1. Therefore, single parti- 
cle distributions of produced secondaries, if only measured 
very precisely, can a priori provide us information not only 
on the temperature T of the hadronizing system but also, 
when investigated using a nonextensive approach, give us 
information (via value of q — 1) on its total heat capac- 
ity C . In this way one can not only check whether some 
(approximate) thermodynamical state is formed in a single 
collision but also what are its thermodynamical properties 
- a very important feature, especially in what concerns the 
existence and type of the possible phase transitions [2]. 




. 10' 



T = 200 MeV 



(b) 



9C 



9°Q 

#19 



a q = 1.05 
O q = 1 .0 



0.0 0.5 1.0 1.5 2.0 

P T [GeV] 



0.0 0.5 1.0 1.5 2.0 

P T [GeV] 

Fig. 6. (a) Normal exponential pr distributions i.e., q — 1) 
for T = 200 MeV (black symbols) and T = 250 MeV open 
symbols) . (b) Typical event from central Pb + Pb at Ebeam = 
3 A-TeV for = 200 MeV for q = 1 (black symbols) expo- 
nential dependence and q = 1.05 (open symbols). (Reprinted 
from Physica A305, G. Wilk and Z. Wlodarczyk, "Applica- 
tion of nonextensive statistics to particle and nuclear physics" , 
227, Copyright (2002), reproduced by permission of Elsevier, 
ht t p : //www . elsevier . com [ ) . 

The next question is: how plausible is such a program? 
The point is, as discussed above, that one performs fits us- 
ing T and q in a Tsallis distributions rather than only T 
in the usual exponential ones. However, the corresponding 
data on are effectively integrated over the longitudinal 
phase space (or, at least a part of it) and are averaged over 
many events. The best thing would be to observe such an 
effect in a single events, then an event-by-event analysis of 
data would be possibly. Fig. [6] shows what we can expect. 
Two scenarios are demonstrated there: (a) T is constant 
in each event but (because, for example, of different ini- 
tial conditions) it fluctuates from event to event and (b) 
T fluctuates in each event around some mean value To. 
We have chosen for comparison a typical event obtained 
in simulations performed for central Pb + Pb collisions 
taking place for beam energy equal Ef, e arn = 3 A-TeV (ex- 
pected shortly in ALICE experiment at LHC). Density of 
particles in the central region (defined by rapidity win- 
dow — 1.5 < y < 1.5) is chosen to be equal to ^ = 6000. 
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In case (a) in each event one expects exponential depen- 
dence with T = T event and possible departure from it 
would occur only after averaging over all events. It would 
reflect fluctuations originating from different initial con- 
ditions for each particular collision. This situation is illus- 
trated in Fig. [6] (a) where px distributions for T = 200 
MeV (black symbols) and T = 250 MeV (open symbols) 
are presented. Such values of T correspond to typical un- 
certainties in T expected at LHC accelerator at CERN. 
Notice that both curves presented here are straight lines. 
In case (b) one should observe departure from the expo- 
nential behavior already on the single event level and it 
should be fully given by q > 1. This reflects a situation 
when, due to some intrinsically dynamical reasons, differ- 
ent parts of a given event can have different temperatures, 
as we have discussed above. In Fig. [5] (b) black symbols 
represent exponential dependence obtained for T = 200 
MeV (the same as in Fig. [6] (a) ), open symbols show the 
power-like dependence as given by (O with the same T 
and with q = 1.05 (notice that the corresponding curve 
bends slightly upward here) . In this typical event we have 
~ 18000 secondaries, i.e., practically the maximal possible 
number. Notice that here points with highest px already 
correspond to single particles. As one can see, experimen- 
tal differentiation between these two scenarios will be very 
difficult, although not totally impossible. On the other 
hand, if successful it would be very rewarding 0. 

The following remarks are worth to be done at this 
point. 

— We are using rather freely the notion of fluctuating 
temperature. The question then arises whether is makes 
sense. Not going into a detailed dispute, we would only 
like to mention at this point that traces of this idea 
can already be found in [491150] F°l . In particular it is 
important when discussing some peculiarities of the 
phase diagrams, which are important when address- 
ing the question of possible phase transitions between 
QCD and normal matter [50]. What we want to do is 
to bring to ones attention the fact that event-by-event 
analysis allows us (at least in principle) to detect fluc- 
tuations of temperature taking place in a given event. 
This is more than an indirect measure of fluctuations 
of T proposed some time ago in [52] or more direct 
fluctuations of T from event to event discussed in [IS] . 

— As the heat capacity C is proportional to the volume, 
C oc V, in our case V would be the volume of the inter- 
action (or hadronization) , it is expected to grow with 
volume and, respectively, q is expected to decrease 
with V. This is indeed the case if one puts together 
the results for e + e~, pp and pp and and AA collisions 
for example, those of [23] for e + e~ collisions, together 

9 It is interesting to realize that for the Planckian gas at 
T = 186 MeV, occupying volume of the order of the volume of 
sulfur nucleus, one gets C = 34.4 per degree of freedom, which 
leads, using Eq. (|20|l . to q = 1.015 obtained for such system 
for the pt dependence of produced secondaries. 
10 For those interested in discussion on the problem of internal 
consistency (or inconsistency) of the notion of fluctuations of 
temperature in thermodynamics we refer to 48 J49II5T] . 



with those of [6] for pp collisions and all results for 
heavy-ion collisions, like [24,20,28] and especially [21] 
where such a trend was found when analyzing heavy- 
ion events with different centrality (i.e., with different 
volumes V). 

— As the parameter q replaces in some sense the action 
of many not yet identified dynamical factors, one ex- 
pects that, with such factors included, q should dimin- 
ish. This is precisely what has been demonstrated in 
[21] analyzing transverse momenta of pions produced 
in RHIC experiments by using a simple minded Tsal- 
lis formula and an accordingly modified Hagedorn [53] 
approach which already contains in it some dynamics 
(based on a special bootstrap hypothesis of resonances 
composed out of resonances itd.). In the second case 
the values of q — 1 found are much smaller, but still 
remain nonzero indicating therefore the existence of 
some residual additional dynamic there. 

— As demonstrated in [21] . using a nonextensive version 
of the statistical model allows us to describe data well 
in the domain previously believed to be governed en- 
tirely be pure jet physics. Deviations (i.e., dominance 
of truly hard collisions) start at px near 10 GeV and 
further. It would mean that the so called mini-jet re- 
gion can probably also be investigated using a nonex- 
tensive approach (what should be, however, checked in 
more detail int the future). 

2.5 The whole phase space 

Already presenting results for the longitudinal phase space 
we encountered multiplicity distributions of produced sec- 
ondaries. They involve the whole of phase space, both its 
longitudinal and transverse components. However, as we 
have already stressed, because (<?l — 1) >> (qx — 1), the 
dominant role of the longitudinal dynamic in establishing 
the actual number of produced secondaries and its fluctu- 
ation from event to event is obvious. We shall now discuss 
this problem in more detail. 

Previous findings could be summarized in the following 
way: knowing the amount of energy W which is going to 
be transferred to the produced secondaries (i.e., knowing 
the inelasticity K of reaction [39] ) and the mean num- 
ber of produced secondaries, (N), and respecting the fact 
that they are essentially distributed in the longitudinal 
phase space only, one arrives, after using the information 
theory approach (cf., [5]), with the usual exponential dis- 
tribution in E = (nix) cosh y. Additional information on 
the fact that produced secondaries are distributed not ac- 
cording to a Poisson distribution but rather according to 
NB distribution characterized by parameter k is enough 
to get g-exponential distribution in E with q = 1 + 1/k. 

Now, it turns out that the opposite is also true, namely, 
as we have shown in [22], the fact that N particles are 
distributed in energy via A^-particle Tsallis distribution 
described by the nonextensivity parameter q allows us to 
show that their number distribution has to be of the NB 
type with k = l/(q — 1). To illustrate this we first start 
with the derivation of Poisson multiplicity distribution 
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and then to compare it with the corresponding derivation 
of the NB distribution [22l. 



2.5.1 Poisson multiplicity distribution 

This distribution arises in a situation where in some pro- 
cess one has N independently produced secondaries with 
energies {JSi,. ..,#}, each distributed according to the sim- 
ple Boltzmann distribution: 



2.5.2 Negative Binomial multiplicity distribution 

This distribution arises when in some process N indepen- 
dent particles with energies {.Ei .,jv} which are now dis- 
tributed according to Tsallis distribution, 



h({E lt ... iN }) = C 



N 



l-(l-g) 



x 



- + 1-N 



/(E()= l.exp(-§ 



(21) 



with normalization constant Cn given by 
1 N 

C N = - w l[[(t-2) q -( i -3)} 



(28) 



(where A = (E)). The corresponding joint probability dis- 
tribution is then given by: 

/({S 1 ,...,jv}) = ^-expf-i^^ i j . (22) 



(g-1) 
X N 



n r[N + 



2-g 
9-1 



r 



(a) 



(29) 



For independent energies {£?j=i,...,iv} the sum E = Ei 
is then distributed according to the following gamma dis- 
tribution, 



9n(E) = 



N-l 



w=w ■ (! )•-*(-! i ^ 



distribuant of which is 

JV-l 



It means that there are some intrinsic (so far unspecified 
but summarily characterized by the parameter q) fluctu- 
ations present in the system under consideration. In this 
case we do not know the characteristic function for the 
Tsallis distribution, however, because we are dealing here 
only with variables {-Ej=i,...,jv} occurring in the form of 
the sum, E — J2iLi E ii one can still sequentially perform 
integrations of the joint probability distribution (|28[) and, 
noting that (as before, cf. eq. ([25]) ) 



G 



N 



;<i-l)! \J) •«*><-!>• < 24 > 



Eq. ([23)1 follows immediately either by using characteristic 
functions or by sequantially performing integration of the 
joint distribution (|22[) and noticing that: 



h N (E) = h N ^(E)- 



E 



N-l (N 

we arrive at formula corresponding to eq. (|23[) . namely 



(30) 



h N (E) 



(n- iy.x N 



E 



N 



9n{E) = g N -!(E) N _ - . 
For energies such that 



N+l 



(25) 



(26) 



,E 



(31) 

+1-N 



i=0 



i=0 



the corresponding multiplicity distribution has a Poisso- 
nian form (notice that E/X= (N)): 



with distribuant given by 

JV-l 

l-J2Hi(E) 

3=1 
E i-1 



H N {E) 
Hi(E) 



where 



P(N) =G N+1 {E)-G N (E) = 



(27) 



l - A -' exp(-aE) = ^-.eM-(N)). 



f[[(i-l) q -(i-3)} 



l-(l-q) 



(32) 



-1-3 



AT! 



N\ 



In other words, whenever we have variables E\ t _^N,N+i,... 
taken from the exponential distribution / (Ei) and when- 
ever these variables satisfy the condition J2i=o E i < E < 
Sj=o E ii then the corresponding multiplicity N has a 
Poissonian distributional. 

11 Actually, this is the method of generating Poisson distri- 
bution in the numerical Monte-Carlo codes. 



As before, for energies E satisfying the condition given 
by eq. (|2"6"|) , the corresponding multiplicity distribution is 
equal to 

P{N)=H N+l (E)-H N (E) (33) 
and is given by the Negative Binomial distribution : 



P(N) = 



(g-1) 
Nl 



N 



i r[N + i + 2 ^ 
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!-(!-<?) 



A 

r(N + k) 
r(N + i)r(k) 



n -N+-r±- 

l — q 



N 



N+k ■ 



(34) 



+ 

where the mean multiplicity and variance are, respectively, 

(35) 



Var(N) = 



l-(l-«) 



= ( N ) + (N) 2 .(q-l). 



This distribution is defined by the parameter k equal to 

1 



(36) 



Notice that for q — ► 1 one has k — > oo and P(N) becomes 
a Poisson distribution, whereas for q — > 2 one has — > 1 
and we are obtaining geometrical distribution 0- 



3 Further developments 



Introducing a kinetic coefficient of conductance v = r\j p 
and denoting k = c p /cv, where cy is specific heat under 
the constant volume for which 1/cy = q—1, we have that 

T ef/ = T + — /(«) = T + (g - 1)— /(u) (41) 

or, because tD = q — 1, one can write this also as 



T 



e/7 



To + (9 " 1)2 ^D /(u) 



(42) 



In [28j the transverse momentum spectra of pions and 
protons and antiprotons produced in the interactions of 
P + P, D + Au and Au + Au at ^/IWn = 200 GeV 
at RHIC-BNL [S3] were analyzed using a nonextensive 
approach. Among other things they found dependencies 
of the nonextensivity parameter q and temperature T on 
the number of participants, N p (i.e., number of nucleons 
taking part in a given AA collision in the production of 
secondaries). From them we have obtained a dependence 
of T on the parameter q which are shown in Figs. [7] and 
[U In all cases we find that f(u) < and that T seems to 
be linearly dependent on q — 1 . 



Let us now proceed a step further in Eq. 
it in the following form, 



T0]1 by writing 



c p p- 



dT 
dt 



a(T'-T)+r)f(u), 



(37) 



with a new term, rjf(u), which presents the effect of a 
possible viscosity (with viscosity coefficient rf) existing in 
the system. The function f(u) contains terms dependent 
on the velocity in the form of 5^ — . Using as before 

OX fc OX i 

T' defined by ([7]) we get an extension of Eq. (fTTj) : 



dT 

~dt 



i 

T= -T 

T 



T, 



rrr 
+ — 

c P P' 



1 

c p p 

f(u) 



(38) 



This equation leads to the Langevin equation resulting in 
fluctuations of the temperature T given in the same form 
of Eq. (TT2"|) as before but with 



1 



/' : 



1 



T, 



—J\ u ) 

Cpp 



T 



eff 



(39) 



In this way previous T = (T) has now been replaced by 
a kind of effective temperature 



T 



eff 



T 







c p p- 



/(«) 



T 



o 



V 



a 



/(«)• 



(40) 



12 Actually the parameter k in NB can be simply expressed 
by the correlation coefficient p for the two-particle energy cor- 
relations, k — (p + 1) / p, see [22] for details. 



0.14-> 




0.07 0.08 0.09 0.10 0.11 
q-1 

Fig. 7. Dependence of temperature T (in GeV) on the pa- 
rameter q for production of negative pions in different reac- 
tions. The solid line shows a linear fit to obtained results: 
T = 0.22 - 1.25(g - 1) (cf., Eq. (gTJ|) and dashed line shows 
the corresponding quadratic fit: T 

El). 



0.17-7.5(q-l) 2 (cf., Eq. 



These results can be compared with old results for 
e + e _ annihilation reactions discussed some time ago in 
terms of g-statistics in |23] . The q dependence of the tem- 
perature parameter T which can be deduced from them is 
shown in Fig. [9] Notice that now the temperature is lower 
and depends only weakly on q. 

Finally, let us discuss results on fluctuations of multi- 
plicity observed in heavy-ion collisions 55J. They exhibit 
non-monotonic changes as function of the number of par- 
ticipants N p [55] . Actually, also changes of (N) show non- 
linear increase, though not so spectacular. Acting in the 
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0.08 



Fig. 8. The same as in Fig.[7]but for the produced antiprotons. 
The linear fit (solid line) is: T = 0.36 — 3.4(g — 1) whereas 
quadratic one (dashed line) is: T = 0.29 — 35(g — l) 2 . 
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Fig. 9. Dependence of temperature T (in GeV) on the parame- 
ter q for the production of pions in e + e" annihilation reactions. 
T = 0.131 - 0.24(g - 1) . 



spirit of our analysis here we can expect that 



Var(N) 
<N> 



<N> 



but now, with T e f f we can show that 



< N >-- 



W 



T, 



eff 



(43) 



(44) 



where T e ff is given by Eq. ([41]) and where W is the full 
accessible energy. We have therefore that 



(N) - n N p 
<N> 



c(q-l). 



(45) 



Here no is the multiplicity in the single nucleon-nuclcon 

VT f(u) 
K T 



collision measured in the region of acceptance, c 
(notice that c is positive because, as was found from Figs. 
0and[H], f(u) < 0). 

In Figs. [10] and [TTJ we show what can be extracted 
from PbPb collision data taken by NA49 experiment [55] . 
As one can observe the data confirm our expectation that 
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Fig. 10. Comparison of Var(N) / (N) versus N p (squares) with 
1 + c(< N > —noNp) vs N p (circles) (here no = 0.642 and 
c = 4.1 . Data are for negatively charged particles from PbPb 
collisions as collected by NA49 experiment [55] . 
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Fig. 11. The same as in Fig. [10] but translated to q— 1 vs -/V p . 
Squares were obtained from Var(N) / (N) vs N p and circles 
from (JV> vs N p . As before n = 0.643 and c = 4.1. 



dependencies of Var(N) / (N) and (N) on the number of 
participants N p are, after introducing the concept of T e ff, 
essentially the same. The value of n = 0.642 is also sen- 
sible, being only a little greater than the multiplicity ob- 
served in pp collisions when calculated using the accep- 
tance of the NA49 experiment. Notice also that value of 
c = 4.1 obtained here for PbPb collisions is not far from 
the value 1.25/0.22 = 5.7 obtained for data from RHIC 
(i.e., for AuAu collisions but at much higher energy) which 
we have obtained in Fig. [7] 

We close this part by noticing that this problem is 
not trivial since none of the known models for multipar- 
ticle production processes describes Var(N)/ (N) vs N p 
observed experimentally |55j . They are described only by 
some specialized models addressing fluctuations, like the 
percolation model 56J, the model assuming inter-particle 
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correlations caused by the combination of strong and elec- 
tromagnetic interactions [57] or the transparency, mixing 
and reflection model [58] . Actually, all those attempts were 
addressing only Var(N)/(N) vs N p but not (N) vs N p . 
From this perspective, results presented above in Figs. [TU1 
and QT] confirm the reasonableness of the idea of T e / / in- 
troduced in this Section. If one uses Var(N) / (N) vs N p 
to obtain q — 1 then it turns out that the same value of 
q — 1 describes also dependence of (N) on N p ; this can 
only be using T e f f and this is because it depends on q — 1. 



4 Remarks and summary 

Let us start with two remarks which are in order here: 

(i) Results which recall directly to Tsallis entropy were 
obtained using the constraint Y^,Pi = 1 an d the for- 
mula J^iPl^i — (A)q for the g-expectation values. On 
the other hand, there exists a formalism, which ex- 
presses both the Tsallis entropy and the expectation 
values using the so-called escort probability distribu- 
tions [59 : Pi — Pi/YliPi- However, as was shown in 
|60j . such an approach is different from the normal 
nonextensive formalism because the Tsallis entropy ex- 
pressed in terms of the escort probability distribu- 
tions has some difficulty with the property of concav- 
ity. From our limited point of view, it seems that there 
is no problem in what concerns practical, phenomeno- 
logical applications of nonextensivity as discussed in 
the present work. Namely, using Pi one gets distribu- 



tions of the type c [1 — (1 — q)x/l 



19/(1-9) 



which is, in 



fact, formally identical with c [1 — (1 — Q)x/ L] 1 ^ 1 ®\ 
provided we identify: Q = 1 + (q — l)/q, L = l/q 
and c = (2 — Q)/L = 1/7. The mean value is now 
(x) = L/(3 - 2Q) = 1/(2 - q) and < Q < 1.5 (to 
be compared with 0.5 < q < 2). Both distributions 
are identical and the problem, of which of them better 
describes data is artificial. 
(ii) One should be aware that there is still an ongoing dis- 
cussion on the meaning of the temperature in nonex- 
tensive systems. However, the small values of the pa- 
rameter q deduced from data in transverse phase space 
(where the connection with thermodynamical approach 
makes sense, as discussed before) allow us to argue 
that, to first approximation, T can be regarded as the 
hadronizing temperature in such a system. One must 
only remember that in general what we study here is 
not so much the state of equilibrium but rather some 
kind of stationary state. For a thorough discussion of 
the temperature of nonextensive systems, see [61] . 

With the above reservations in mind, we can summa- 
rize that, when looking from the point of view of a sta- 
tistical approach [HI], the power-law behavior of many 
distributions observed in elementary and heavy ion col- 
lisions can be traced back to the necessity of using the 
nonextensive version of a statistical model (here taken in 
the form proposed by Tsallis [5]). 

We interpret this as a sign of some intrinsic fluctua- 
tions present in any hadronizing system, which were only 



recently to be recognized as vital observable when search- 
ing for the production of the QGP form of matter [2]. In 
fact, a number of works [62j have demonstrated the ex- 
istence in such reactions of event-by-event fluctuations of 
the average transverse momenta (p) per event. The quanti- 
ties considered were: Var ((p)) / ((p)) 2 and (ApiApj)/ ({p}} 2 
These quantities can be shown 22J to be fully determined 
by uj as defined by Eq. (fj|, i.e., by fluctuations of the 
temperature T of hadronizing system - a vital observable 
when searching for QGFP^I. In fact, when considering the 
case of N ev events with Nk particles in the fc th event, one 
has that 

Var((p)) _Var(T) 



(TY 



where 



«P» = 



N K 



(p) k ; with (p) fe = — Pi > (47) 



This is what we have shown in the last part od Section [3] 
This is the problem which needs further investigations. 
We close with some remarks: 

— Although our original investigations presented here were 
based on the notion of Tsallis entropy (usually with the 
help of information theory) one must mention that one 
can also get Tsallis distribution without resorting to a 
Tsallis entropy altogether (see, for example, [64]). 

— The other way to get a Tsallis distribution from some 
general thermodynamical considerations was presented 
in [65j . It is based on allowing a linear dependence 
of the temperature T on energy, T = T + (q — 1)E. 
Here temperature is not fluctuating. Actually, if one 
would like to follow this approach and to have Tsal- 
lis distribution with T e f f discussed in Section [3] one 
should write T = T + (q - 1) • const + (q - 1)E. 
Then T = Tq only would result in exp (— E/Tq), T = 
To(q — V)E would result in the usual Tsallis distribu- 
tion cxp q {-E/Tq), T = T eff = T + (q - 1) • const 
would give exp (— E/T e ff] and, finally, T — T e ff — 
T +(q-l)-const+(q-l)E would giveexp^ (-E/T eff ). 

— Notice that for x >> X/(q— 1) Tsallis distribution 
becomes a pure power low and loses its dependence on 
the scale A: f(x) - [l-(l-q)x/X] 1 ^ 1 -^ — ► ^/(l-s). 

— Instead of using an intrinsic fluctuations one can also 
obtain a power law distribution by using the notion 
of self-organized criticality [SH] in a cascade processes 

( C f., n5Ki). 

— Another interesting possibility, not yet fully explored, 
is that, as shown in [68] , one can formulate a descrip- 
tion of the so called stochastic networks using nonex- 
tensive information theory based on Tsallis statistics. 
Using this approach one can then demonstrate [61?] 



13 Generally speaking, an analysis of transverse momenta pr 
alone indicates very small fluctuations of T. On the other hand, 
as reported in [63] , the measured fluctuations of multiplicities 
of produced secondaries are large (i.e., multiplicity distribu- 
tions are substantially broader than Poissonian). 
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that hadron production viewed as formation of a spe- 
cific stochastic network can explain in a natural way 
the power-law distributions of transverse mass spectra 
of pions found in [70] . 

— In the string models of production of hadrons the nat- 
ural distribution in pr is exp (— 7rmy/re) rather than 
exp (mr/T) (where re is string tension) really observed. 
However, if one allows for the gaussian fluctuations 
of the parameter re (characterized by parameter (re 2 ) 
which can be connected with the fluctuations in the 
QCD vacuum) then the first form is transformed into 
the second one with T = y/ (re 2 ) / (2tt) (i.e., in this ap- 
proach parameter T characterizes rather the properties 
of the QCD vacuum than those of hadrons) - see [71] . 

— Finally, recall that when applied to the hydrodynam- 
ical model of multiparticle production the nonexten- 
sivity approach converts the usual nonviscous hadronic 
fluid into the viscous one preserving, however, the usual 
linear flow equations (albeit now given in a nonexten- 
sive form) [72]. 

Partial support (GW) of the Ministry of Science and Higher 
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